集合卡尔曼滤波同化中雷达位置的敏感性研究

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。      

EnSRF雷达资料同化在一次飑线过程中的应用研究

本文利用包含复杂冰相微物理过程的WRF（Weather Research and Forecasting）模式,针对2007年4月23日发生在我国华南地区的一次典型飑线天气过程,分别进行了确定性预报和集合预报试验,发现确定性预报能大致捕捉到飑线系统的发生发展过程,但对飑线后部的层云区模拟效果较差。集合预报能够有效地减少模式的不确定性,大部分集合成员对飑线的模拟效果优于确定性预报。进一步将集合预报得到的40个成员作为背景场,采用EnSRF（Ensemble Square Root Filter）同化多普勒天气雷达资料,并将分析得到的集合作为初始场进行集合预报,通过与未同化雷达资料的集合对比,考察了EnSRF同化多部雷达资料对飑线系统的影响。结果表明:EnSRF雷达资料同化增加了模式初始场的中小尺度信息,大部分集合成员的分析场能够较准确地再现飑线的热力场、动力场和微物理场的细致特征,并且模拟出飑线后部的层云结构。通过对EnSRF分析的集合进行模拟发现,大部分集合成员较未同化雷达资料时模拟效果有明显改善。同化后的集合预报ETS（Equitable Threat Score）评分最高,其次是未同化的集合预报,确定性预报的最低。     

中低纬太平洋夏季次表层海温异常年代际变化的时空特征比较

对1955—1998年的夏季次表层(0~400 m)海温进行了EOF分解,对比分析了中低纬太平洋夏季次表层海温的年代际变化特征。结果表明:中纬太平洋夏季次表层海温有2种年代际分布模态:0~160 m表现为PDO型,200~400 m表现为全区一致型;低纬太平洋夏季次表层海温有3种年代际分布模态:0~60 m和300~400 m为全区一致型,80~240 m为东西反向型。太平洋夏季次表层海温异常的年代际变化在中低纬都存在从上而下的时间滞后;而同一层中低纬太平洋夏季次表层海温年代际突变的时间也不一致。     

Hybrid ETKF-3DVAR方法同化多普勒雷达速度观测资料Ⅰ:模拟资料试验

利用WRF(Weather research and forecasting)模式及模式模拟的资料,采用Hybrid ETKF-3DVAR(ensemble transform Kalman filter-three-dimensional variational data assimilation)方法同化模拟雷达观测资料。该混合同化方法将集合转换卡尔曼滤波(ensemble transform Kalman filter)得到的集合样本扰动通过转换矩阵直接作用到背景场上,利用顺序滤波的思想得到分析扰动场;然后通过增加额外控制变量的方式把"流依赖"的集合协方差信息引入到变分目标函数中去,在3DVAR框架基础下与观测数据进行融合,从而给出分析场的最优估计。试验结果表明,Hybrid ETKF-3DVAR同化方法相比传统3DVAR可以提供更为准确的分析场,Hybrid方法雷达资料初始化模拟的台风涡旋结构与位置比3DVAR更加接近"真实场",对台风路径预报也有明显改进。通过对比Hybrid S试验与Hybrid F试验发现,Hybrid的正效果主要来源于混合背景误差协方差中的"流依赖"信息,集合平均场代替确定性背景场带来的效果并不显著。     

历史数据和强化学习相结合的低频轨迹数据匹配算法

针对低频(采样间隔大于1min)轨迹数据匹配算法精度不高的问题,提出了一种基于强化学习和历史轨迹的匹配算法HMDP-Q,首先通过增量匹配算法提取历史路径作为历史参考经验库;根据历史参考经验库、最短路径和可达性筛选候选路径集;再将地图匹配过程建模成马尔科夫决策过程,利用轨迹点偏离道路距离和历史轨迹构建回报函数;然后借助强化学习算法求解马尔科夫决策过程的最大回报值,即轨迹与道路的最优匹配结果;最后应用某市浮动车轨迹数据进行试验。结果表明:本文算法能有效提高轨迹数据与道路匹配精度;本算法在1min低频采样间隔下轨迹匹配准确率达到了89.2%;采样频率为16min时,该算法匹配精度也能达到61.4%;与IVVM算法相比,HMDP-Q算法匹配精度和求解效率均优于IVVM算法,16min采样频率时本文算法轨迹匹配精度提高了26%。     

农业遥感研究应用进展与展望

得益于中国自主遥感卫星、无人机遥感和物联网等技术的发展,中国农业遥感研究与应用在过去20年取得了显著进步,中国农业遥感信息获取呈现出天地网一体化的趋势;农业定量遥感在关键参数遥感反演技术方法与应用方面取得进展;作物面积、长势、产量、灾害遥感监测的理论与技术方法取得突破,农业遥感技术应用领域不断拓展。本文从农业遥感信息获取、农业定量遥感、农业灾害遥感、作物遥感识别与制图、作物长势遥感监测与产量预测、农业土地资源遥感等方面对中国农业遥感科研与应用进行了总结综述。     

胶东玲珑金矿田175号脉群深部地质特征与成矿预测

玲珑金矿属于以石英脉型为主的金矿床,其矿脉受玲珑断裂和破头青断裂控制。区内脉群发育,175号脉群是玲珑金矿的主要大型金矿脉群之一。文章通过总结175号脉群成矿地质背景、矿床地质特征,分析了控矿因素、矿体赋存规律、矿石类型等;阐述了175号脉群深部成矿规律,并根据见矿工程利用地质块段法预测175号脉群11—101勘探线范围-800 m标高以下(332+333)金矿石量192.44×10~4 t,金金属量8.5t。     

一幅新的第四纪地质及地貌图的编制——以宁夏与山东为例

为了更好地与生态环境及灾害地质紧密结合,第二代全国区域地质志编制项目规定第四纪地质及地貌图为必编图件之一。图件内容包括第四纪地质体的时代、成因类型、岩性和地貌成因形态组合类型、活动断裂及海侵范围等重要地质内容,同时,附第四纪地质-地貌典型剖面和重要的地貌景点。图件特点:(1)图件的地理底图是首次用数字高程模型(DEM)数据灰度图作背景;(2)第四纪地质内容与地貌成因形态类型同编为一幅图;(3)第四纪地质体的面色用成因类型色表示;(4)地貌部分划分了不同级别成因形态类型和有特色的微地貌景点。总之,图面有立体感,地形高差明显,层次分明,色彩鲜艳,是一幅崭新的图件,可供防灾治理和地质旅游参考应用。     

我国海岛管理信息化建设问题与对策研究

海岛管理信息化建设是《中华人民共和国海岛保护法》赋予海洋行政主管部门的主要职责之一。文章分析研究我国海岛管理信息化建设的背景和现状,指出面临的一系列问题,从总体规划、数据应用、完善机制、统筹管理、加强配套、保障安全等方面提出相应对策,为我国海岛管理信息化建设提供参考。     

The effects of clay minerals and organic matter on nanoscale pores in Lower Paleozoic shale gas reservoirs,Guizhou, China

In organic-rich gas shales, clay minerals and organic matter (OM) have significant influences on the origin, preservation, and production of shale gas. Because of the substantial role of nanoscale pores in the generation, storage, and seepage of shale gas, we examined the effects of clay minerals and OM on nanoscale pore distribution characteristics in Lower Paleozoic shale gas reservoirs. Using the Niutitang and Longmaxi shales as examples, we determined the effects of clay minerals and OM on pores through sedimentation experiments. Field emission–scanning electron microscopy combined with low-pressure N₂ adsorption of the samples before and after sedimentation showed significant differences in pore location and pore size distribution between the Niutitang and Longmaxi shales. Nanoscale pores mostly existed in OM in the Longmaxi shale and in clay minerals or OM–clay composites in the Niutitang shale. The distribution differences were attributed largely to variability in thermal evolution and tectonic development and might account for the difference in gas-bearing capacity between the Niutitang and Longmaxi reservoirs. In the nanoscale range, mesopores accounted for 61–76% of total nanoscale pore volume. Considerably developed nanoscale pores in OM were distributed in a broad size range in the Longmaxi shale, which led to good pore connectivity and gas production. Numerous narrow pores (i.e., pores?<?20 nm) in OM–clay composites were found in the Niutitang shale, and might account for this shale’s poor pore connectivity and low gas production efficiency. Enhancing the connectivity of the mesopores (especially pores?<?20 nm and those developed in OM–clay composites) might be the key to improving development of the Niutitang shale. The findings provide new insight into the formation and evolutionary mechanism of nanoscale pores developed in OM and clay minerals.